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ANALYTIC SOLUTIONS TO THE IGNITION KINETICS 


OF THE HYDROGEN-OXYGEN REACTION 


by Richard S. Brokaw 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

SUMMARY 

Solutions have been obtained assuming isothermal conditions and neg- 
ligible depletion of reactants during the induction period for the follow- 
ing scheme of reactions: 

k l 

OH + H 2 H 2 0 + H 

k 2 

H + 0 2 ►OH + 0 



0 + 



OH + H 


k 6 

H + 0 2 + M ► H0 2 + M 

k ll 

H0 2 + H 2 * H 2 0 2 + H 


Three types of solutions have been deduced, corresponding to kgc^ > 2k 2 
(a low-temperature - high-pressure region of long ignition delays), 

2k 2 > kgcM (a high-temperature - low-pressure region of short delays), 
and 2k 2 = kgCj^ (the boundary between the two regions). Experimental 
ignition delays in both long and short delay regions are adequately ex- 
plained, and ignition delays computed by numerical integration of the 
rate equations are reproduced to with a few percent. Finally, ignition 
delays for the Von Neumann spike condition in Chapman- Jouguet detonations 
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in the neighborhood of the learn limit of detonability in air have been 
calculated. These delays decrease by two orders of magnitude as the hy- 
drogen mole fraction increases from 0.135 to 0.145. This is very close 
to the experimental limit (hydrogen mole fraction 0.14 to 0.15) and thus 
supports the notion that the condition 2kg > kgCj^ must be satisfied 
for a stable detonation. 

INTRODUCTION 


When hydrogen and air are instantaneously mixed at high temperature 
(as in idealized analyses pertinent to hypersonic ramjet s^> ^), or when a 
mixture containing hydrogen and oxygen is subjected to a pulse of high 
temperature and pressure (as in shock tube ignition studies > ), the fol- 
lowing events ensue: First, there is a very short period during which 

small concentrations of atoms and radicals - H, OH, and 0 - are produced. 
The mechanism of this initiation is uncertain and may vary under differ- 
ent circumstances. There follows the bulk of the induction period during 
which the chain reaction is propagated and branched through the elementary 
reactions 


OH + H 2 
H + O2 

0 + h 2 


k l 


+-H 2 0 + H 
* OH + 0 

+ 0H + H 


I 

TT 

III 


The first reaction is not strongly exothermic (i.e., 15 kilocalories at 
1000° K), while the second and third reactions are endothermic by 16 and 2 
kilocalories, respectively. Thus, during most of the induction period, 
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the atom and radical concentrations increase exponentially under essen- 
tially isothermal conditions. Very late in the induction period, atom 
and radical concentrations become so large that self-heating occurs, 
largely due to the highly exothermic three-body recombination processes. 
The temperature rises precipitously, and the induction period is at an 
end. 

At the lower temperatures and higher pressures the further reactions 

*6 

H + 0 2 + M * H0 2 + M VI 


H0 2 + H 2 - 


‘n 


■*H 2 0 2 + H 


XI 


must be considered as well. These five reactions, together with terms 
describing the destruction of chain carriers after diffusion to the walls, 
suffice to explain qualitatively the ignition behavior of the hydrogen- 
oxygen system, including the existence of the first, second, and third ex- 
plosion limits. 

The differential equations describing the chemical kinetics of reac- 
tions I, II, and III have been integrated numerically for hypersonic ram- 
.,12 4 

jet > and shock tube conditions. This paper presents analytic solutions 
that are excellent approximations to the numerical results; furthermore, 
reactions VI and XI have been included. These analytic solutions are 
valuable in that they permit one to determine easily the effects of varia- 
tions in specific rate constants, temperature, pressure, initiation 
mechanism, etc., on ignition delays and the like. Some illustrative ap- 
plications are presented and discussed. 
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THEORETICAL CONSIDERATIONS 


The differential equations governing the growth 
tions during the induction period are as follows: 


dc, 


OH 


dt 


= " k l c H 2 c OH + k 2 c 0 2 c H + k 3 c H 2 c 0 


of radical concentra- 


+ i 0H (1) 


dc H 

dt 


= kic H2 c 0H - k 2 c 02 c H + k 3 c H2 c 0 - k 6 c 02 c M c H + k n c H2 c H o 2 + 




( 2 ) 


dc 0 

— = k 2 c 02 c H - k 3 c H2 c 0 + i 0 


(3) 


dc H02 

- = k 6 c 02 c M c H - k i;L c H2 c H o 2 + iho 2 


(4) 


Here, the k's are the specific reaction rate constants, the c’s are 
concentrations, and the i’s are rates at which the radicals are generated 
spontaneously, either in the gas phase or on the surfaces. (It is also 
possible to include destruction of radicals at the walls by including 
terms like Yw c OH> e ^ c ’ 4n Eqs. (l) to (4), with the coefficient Yv a 
function of diffusivity, pressure, geometry, and surface activity. ) 

During most of the induction period, the kinetics of the elementary 
reactions I to III, VI, and XI can be represented as pseudo first-order 
reactions; that is, they are all first order in the atom and radical con- 
centrations - the only concentrations that are changing appreciably. 

The initiation rates iQjj, etc., which are zero order in the radical 
concentrations, can be eliminated from Eqs. (l) to (4) by means of the 
substitutions 5 


C 0H 


i H + i 0 + %) 2 

C ° H+ ^ 


(5) 


5 - 


C H = C H 


2k 2 c 0 2 


, iOH + iH + iH0 2 

c O “ c O + * 


2k 3 c H . 


( 6 ) 

(7) 


c ho 2 = c H0 2 + 


Vm 


2k 2 k ll c H 2 


] 


iQH + % + 2i O + 



These initiation rates are presumably very small compared with the rates 
of reactions I to III, so that the correction terms of Eqs. (5) to (8) are 
negligible except in the very initial portion of the induction period. 
Consequently, in what follows, these corrections will be neglected. They 
should be considered, however, in any detailed analysis of the kinetics 
of initiation. Inclusion of these terms poses no mathematical problem - 
they simply increase algebraic complexity. 

The system of differential equations (Eqs. (1) to (4)) may be solved 
by the method outlined by Frost and Pearson®, assuming that the concentra- 
tions of reactants - H 2 , 0 2 , and M - are constant and neglecting iQg, ig, 
io, and iH0 2 ( or applying Eqs. (5) to (8)). A set of particular solutions 
is 

c ± = A i e At i = OH, H, 0, H0 2 (9) 

Substituting Eqs. (9) into Eqs. (l) to (4) yields the relations 



~( k l c H? + 


+ k 2 c 0 2 A H + k 3 c H 2 A 0 


= 0 
( 10 )' 


k l c H 2 A 0H “ ( k 2 c 0 2 + k 6 c 0 2 c M + 


+ k 3 c H 2 A 0 


+ kn t C 


11 c H 2 A H0 2 ■ 


0 

( 11 ) 


k 2 c 0 2 A H " ( k 3 c H 2 + ^^0 


= 0 
( 12 ) 


k 6 c 0o c M A H 


- ( k ll c Hp + A ) A H0 2 = 0 


(13) 


A nontrivial solution requires that the determinant of the coefficients of the 
he zero®. This leads to the quartic equation 

A 4 + [(k-L + k 3 + k l;L )c H2 + (k 2 + k 6 c M )c 02 jA 3 

+ [ k l k 3 c H 2 + ( k l + k 3^ k 6 c H 2 c 0 2 c M + ( k l c H 2 + k 2 c 0 2 + k 3 c H 2 ) k ll%] a2 

+ k i k 3 c H 2 |( k 6 c M " 2k 2^ c 0 2 + k ll c H J A " 2k l k 2 k 3 k ll c H 2 c 0 2 = 0 ( 14 ) 

We seek the positive roots of this equation, which govern the exponential growth 
of radical concentrations during the induction period. The negative roots are 
transients important only in the initiation phases of the reaction - they 
quickly die out. 

The rate constant k -, is very small in comparison with k-^, kg, k^, and 
k 6 c M* Therefore, the additive terms involving k-Q in the coefficients of A, 

O 3 

A , and A may he neglected. 


We shall derive positive roots of Eq. (14) for three conditions: 
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k 6?M > 2k 2* This re S ime corresponds to the conditions above the third ex- 
plosion limit. This is a region of long ignition lags, with delays of 
several seconds at 560° to 600° C and atmospheric pressure. Correspond- 
ingly, the positive A is very small and we find that A 2 , A 3 , and A 4 
may he set equal to zero in Eq. (14). Thus, we obtain 

_ 2k 2 k ll c H 2 

A = k 6 c M - 2k 2 (15) 


This expression is a very good approximation indeed, except very close to 
the second limit condition where the denominator approaches zero. Note 
that, at high pressure, A -*■ 2k 2 kjj_cu /kgCj^, which is independent of pres- 
sure. 

k 6 c M = 2k 2 (Second limit condition). Above the junction of the second and 

third limits, this condition marks the boundary between regions of long 

and short ignition delay. Here, the coefficient of A in Eq. (14) is 

zero (neglecting the term in k -[ -[ ). We may neglect A^ and A 4 rela- 
2 

tive to A to obtain 


A = 


2k l k 2 k 3 k ll c H 2 c 0 2 


“I 1 / 2 


k l k 3 c H 2 + ( k l + k 3 ) k 6 c 0 2 c M 


(16) 


2k 2 > kgCpj. This corresponds to the important region of short ignition 
delays at high temperatures and modest pressures. Under these conditions, 
the coefficient of A in Eq. (14) becomes negative, all terms involving 
k xl may be neglected, and Eq. (14) reduces to the cubic equation 

ki k 3 (2k 2 - k 6 c M )c| c 0 

A = 5 — 2 _ 2 

k l k 3 c H 2 + ( k l + k 3) k 6 c H 2 c 0 2 c M + L( k l + k 3) c H 2 + ( k 2 + k 6 c M) c 0^ + A 


( 17 ) 



8 


Note that, if A is small, the terms involving A and A 2 in the de- 
nominator are unimportant. If, in addition, k]_, k 2 , and k 3 are much 
larger than kgC M , A approaches 2k 2 CQ . This provides the basis of 
the ignition delay correlation used by Schott and Kinsey. 

Equation (17) has a further interesting aspect. Suppose we have an 
approximate value of A, which we designate A]_. We can use the right 
side of Eq. (17) to generate a second approximation, A 2 . We seek a 
positive root. Hence, if A^_ is too small, A 2 will be too large; if 
Aj_ is too large, A 2 is too small. Thus, the correct solution is 
bracketed, and the approximation A a (AiV 1 / 2 should be in error by 
less than half the error of the worst of the two estimates. 

An approximate solution can be obtained by invoking the steady- 
state approximation in OH (setting Eq. (l) equal to zero). Thus, 

Eq. (1) can be used to eliminate the OH concentration from Eqs. (2) 
and (3), which can then be solved with the result 


k,CTi + k-Cr. c. 


An = 


V 3 U H, 


6 0 2 M 


1 + 


4k 3 (2k 2 - kgC^Cj^CQ^ 
( k 3°H 2 + k 6 c 0 2 c M^ 2 


1/2 


(18) 


This, in turn, may be substituted into Eq. (17) to obtain A 2 , so that 
the final approximation to A is: 


*Note that if the error in A 2 is less than twice the error in Aj_, 
the geometric mean value (call it A 3 ) is clearly superior to A]_. Thus, 

A 3 can be substituted into Eq. (17) to obtain an improved approximation, 
A^; A 5 = (A^) 1 / 2 will be an improvement over A 3 - (Ai_A 2 )'*'/ 2 . This 
provides an iterative procedure that will converge on the true value of A. 
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It is possible, of course, to obtain other first approximations Ap by as- 
suming steady states in H or 0 rather than OH. Some numerical comparisons with 
these other techniques are shown in table I (taking k g = 0). The assumption of 
a steady state in OH gives clearly superior values of Ap although the final 
approximations ( A]_ Ag ) 1 / 2 are remarkably good in all cases. The OH concentra- 
tion is smallest, so it seems reasonable that the steady- state approximation 
should apply best to this radical. (The true values of A were obtained by 
numerical solution of Eq. (17).) 

Once the value of A has been obtained, relations among Aqjj, Ajj, Aq, and 
A ho „ follow from Eqs. (10) to (13). 

£ j 

It remains merely to discuss the boundary conditions that define the 
beginning and end of the ignition delay period. We shall assume for simplicity 
that the initiation phase of the induction period is very short and may be 
neglected in comparison with the total span of the delay before ignition. (As 
indicated heretofore, it is possible to obtain analytic solutions without this 
assumption, but at the expense of algebraic complexity. Further, a detailed 
knowledge of the initiation kinetics would be required. ) 


« 

% 
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Let us assume that the initiation process produces radicals by some 
bimolecular process: 



( 20 ) 


•where c M represents the total concentration of stable species. Examples 
of reactions of this sort could include 

H 2 + 0 2 *20H 

H 2 + 0 2 — >-H + HO 2 
H 2 + M *2H + M 


and others as well. It is reasonable to assume that, initially, the rate 
of formation of radicals is equal to the rate of consumption in reactions 
such as I, II, and III, etc., (which are first order in the radical, but 
second order overall): 



where is an effective initial radical concentration (see (Eq. 

Thus, in general, from Eqs. (20) and (21), 

Aj_ = const c M = a-^p 


( 21 ) 

(9)). 


( 22 ) 


so that the radical concentration as a function of time is 

Ci = a iP e At (23) 

as long as initiation is in the gas phase and requires only a small frac- 
tion of the ignition delay. 

The end of the induction period can be defined in a number of ways 

•z 

that are qualitatively equivalent. For example, Schott and Kinsey 0 define 


11 - 


the end of the delay period as the time at which the concentration of some 
species reaches a characteristic level (in their case, c QH = 10" 6 moles 
liter" 1 ). Thus, from Eq. (23), 

At = B - In p (24) 

where t is the ignition delay and B is a constant. The same result 
is obtained if one characterizes the end of the delay period by a critical 
rate of temperature rise due to three-body radical recombination: 

C M°p(f ) cr . t * c i c j c M M ij * c MP 2e2AT ( 25 ) 

which again leads to Eq. (24). Again, if one assumes a critical rate of 
temperature rise due to a bimolecular process such as reaction I, Eq. (24) 
follows once more. Thus, a variety of reasonable assumptions about the 
initiation and termination of the induction period all lead to expressions 
of the form of Eq. (24). 

ILLUSTRATIVE APPLICATIONS 

In this section, the foregoing theoretical results are applied to 
calculate ignition delays for hydrogen-air systems under a variety of 
conditions of practical interest including (l) the low- temperature region 
of moderate pressures, where ignition delays of one to several seconds 
are observed, (2) the higher- temperature low-to-moderate-pressure regime, 
where delays range from several microseconds to several milliseconds, and 
finally (3) the high-temperature - high-pressure boundary between these 
regimes, which seems to be of importance in determining limits of 
detonability. 
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Region of long delays . Ignition delays for hydrogen in air -were measured _ 
some three decades ago by Dixon using an apparatus in which a heated 
stream of hydrogen was injected into a larger concentric stream of heated 
air. Interpolated delays at 600° C are shown in Fig. 1. 

In this regime, the characteristic parameter A (used in Eq. (24) to 
compute the ignition delay t) is given by Eq. (15). From Fig. 1 it is 
clear that, at a pressure close to one-half atmosphere, the ignition de- 
lays become very short, corresponding to the condition kgCj^ = 2kg. This 
provides a key for determining the effective hydrogen concentration (which 
is otherwise indeterminate in Dixon's experiment), because kg is a func- 
tion of the nature of the third body, M. We assume 

kg = 3. 27x10 15 (xjj 2 + 0.35 Xq 2 + 0.43 X [| 2 +0.2 x^ 

+ 14.3 Xjj 2 q + . . .Jt " 1 * 92 liters 2 mole" 2 sec -1 (26) 

where the x's axe the mole fractions of various third bodies. Equa- 
tion (26) is consistent with Clyne's® room temperature measurement (with 

Q 

argon as a third body), and Baldwin^ estimates at 540 C. The relative 
third-body efficiencies are taken from reference 9. Using k 0 from 
Eq. (26) and kg from reference 4 we find that the effective mole frac- 
tion of hydrogen is 0.22. This procedure may seem somewhat arbitrary, but 
it is comforting that the hydrogen mole fraction did not turn out to be 
negative (limit pressure > 0.66 atm) or greater than unity (limit 
pressure <0.27 atm). T 

O -| 

The curve in Fig. 1, computed assuming k-^ = 5. 2x10 liters mole 
sec - ^ and B =3.8, provides a very good fit of the experimental data. This 
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value of k-^ is greater than that suggested by Voevodsky^® by a factor 
of 4.5. The cause of this discrepancy is uncertain: Voevodsky perhaps 

deduced his k-^ assuming somewhat different values of kg and kg; but, 
on the other hand, our kinetic scheme is probably oversimplified. 

Other observations by Dixon” 7 are also easily understood from these 
theoretical relationships. Thus, he reported ignition delays in oxygen 
slightly shorter than those in air. This is not surprising since the 
oxygen concentration does not appear explicitly in Eq. (15); its only 
effect is to cause a slight decrease in kg (see Eq. (26)). On the other 
hand, Dixon reported markedly longer delays in the presence of water 
vapor. This is due to the fact that water is a very efficient third 
body, so that kg is considerably increased. 

The effect of hydrogen concentration on the ignition delay has been 

determined by Anagnostou, Brokaw, and Butler; ^ some of their results 

near 600° C are shown in Fig. 2. (They found some variation according 

to the mode of operation of their apparatus; the shortest delays at this 

temperature were roughly half as large as those shown in Fig. 2. ) The 

curve in Fig. 2 has been calculated from Eqs. (15) and (24) with the 

same constants used to fit Dixon's data in Fig. 1. Thus, the experimental 

_o 0 

composition dependence is rather well predicted: t « p^ * . 

2 

It seems then that the analytic solutions presented in this paper 
describe the behavior of the hydrogen-oxygen reaction in the region of 
long delays rather well. 

Region of short delays . Ignition delays for hydrogen- air have been calcu- 
lated numerically by Momtchiloff, Taback, and Buswell, ^ and Belles and 
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Lauver 4 using different assumptions as to the mechanism of initiation. 
Momtchiloff, Taback, and Buswell use Schott and Kinsey's 3 criterion for 
ignition (cqjj = 10" 6 moles liter"-*-); Belles and Lauver 4 give constants 
for Eq. (9) so that Schott and Kinsey's criterion is easily applied to 
their calculations as well. (Both papers neglect reaction VI. ) 

The analysis of reference 1, presumably pertinent to hypersonic 
ramjets, assumes equilibrium concentrations of hydrogen and oxygen atoms 
in the hydrogen and air streams prior to instantaneous mixing. The hy- 
drogen atoms are expected to be more important because of their higher 
equilibrium concentration. Thus, the growth of OH radical concentration 
with time is 


C 0H = 


Adh 

c H,Eq a h 


.At 


(27) 


From Eqs. (10 ) and (ll), neglecting reactions VI and XI, we obtain 

%I 2k 2 c 0 2 * A 

a H 2k l c H 2 + A 


(28) 


Thus, the induction times computed numerically in reference 1 can be ap- 
proximated from Eqs. (19), (27), and (28). Numerical and approximate 
analytic results are compared in table II. 

Belles and Lauver, 4 on the other hand, assume initiation via the 
reaction 

koH 

H 2 + o 2 *20H 

Assuming that, initially, the rate of production of OH radicals in this 
reaction equals the rate of their consumption in reaction I, we obtain 

2k 0H c 0 2 


( 29 ) 


so that the ignition delay is 



Numerical 4 and analytic results from Eq. (30) are compared in ta- 
ble II. It is clear that the analytic solutions are very good approxima- 
tions to the numerical results, both for very lean and stoichiometric 
conditions, and for two rather different assumptions as to the initial 
source of radicals. 

By studying stationary hydrogen- air flames in a flow system, Foure^ 2 
has obtained information on ignition delays in the region of short delays. 
Foure heated flowing homogeneous fuel- air mixtures to a temperature some- 
what below the point of spontaneous ignition. The mixture, traveling at 
a high subsonic velocity, then passed through a diffuser into a reaction 
tube. The velocity decrease to about 40 meters per second was accom- 
panied by increases in temperature and pressure that rendered the mix- 
tures spontaneously inflammable. Temperature and composition profiles 
were obtained by probing the extended reaction zone formed some distance 
downstream of the diffuser. The distance from the diffuser exit to the 
point of appreciable temperature rise provides an indication of ignition 
delay times; however, the choice of the -diffuser exit as a zero point is 
somewhat arbitrary since reaction may actually start within the diffuser. 
Foure' s results are shown in Fig. 3. A millisecond scale has been added 
based on the assumption that the velocity in reaction tube was 40 meters 
per second. 

Ignition delays computed from Eqs. (19) and (24) are shown by the 
curve in Fig. 3. The constant B has been taken as 25 so as to fit the 
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experimental data at the leaner compositions. The deviation of the two 
richest mole fractions suggests that there may have "been considerable 
reaction in the diffuser under these conditions. 

Foure was unable to observe any temperature rise in his ignition 
tube below 910° K. Since the tube was 1 meter long and the flow velocity 
was about 40 meters per second, this suggests an ignition delay in ex- 
cess of 25 milliseconds. Assuming again that B = 25, we compute an 
ignition delay of 54 milliseconds. At this condition the delay is ex- 
tremely sensitive to the temperature - if the temperature is raised but 
1°, the delay decreases to 40 milliseconds, and we would hence predict 
a 25-millisecond delay at 912° or 913° K. Thus, the present analysis 
predicts Foure' s lower temperature limit within 3° K. 

One further point merits some discussion. The ignition delays in 
the region of long delays were fit assuming B 2 3.8 at 600° K, whereas 
the short delay data require B ss 25 at a temperature some 60° higher. 
From Eqs. (20) to (25) it appears that 

/ dc i\ 

B = -In -tt— ) + 2 In p + const 

\ dt /t=0 

so that a small value of B indicates a large initiation rate and vice 
versa. Thus, it appears that the initiation mechanisms in the long and 
short delay regions are distinctly different, and the long delay initia- 
tion rate is larger by perhaps nine orders of magnitude (although this 
changes the delays by less than one order of magnitude). There is evi- 
dence^ that hydrogen peroxide plays an important role in the initiation 
of the hydrogen-oxygen reaction in the neighborhood of the second limit. 
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Hydrogen peroxide is produced in reaction XI, which is important prin- 
cipally in the long delay regime. 

Boundary between long and short delay regions (and limits of detonability) . 
A stable detonation wave involves the coupling of a shock front with a 
rapid chemical heat release, which, in turn, sustains the shock. Thus, 
the heat release must occur rather close to the shock in a stable 
detonation. 

Belles-*- 4 has proposed that, in the case of the hydrogen- oxygen reac- 
tion, the limits of detonability may be predicted by requiring that, in 
a stable detonation, the second limit condition 2k 2 > kgCj^ be satisfied 
in the Von Neumann spike accompanying the detonation front. (These con- 

ditions of temperature and pressure are calculable from the Chapman- 
Jouguet detonation Mach number, and the adiabatic shock relations for the 
nonreacting gas at the same Mach number. ) Using the second- limit criterion, 
Belles predicted lean and rich limits of detonability for hydrogen- oxygen 
and hydrogen-air in good agreement with experiment. 

Belles applied the second-limit criterion under conditions of tem- 
perature (~1300° K) and pressure (~20 atm) where hydrogen- oxygen and 
hydrogen- air mixtures invariably ignite spontaneously (in other words, at 
temperatures above the junction of the second and third explosion limits). 
Thus, he tacitly assumes that there is a "memory" of the second limit - a 
demarcation between regions of long and short delay. 

To examine this notion further, ignition delays have been calculated 
for conditions around the lean limit of detonability for hydrogen-air mix- 
tures. Chapman- Jouguet Mach numbers were computed as described by Zeleznik 
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and Gordon, 15 and then conditions immediately behind the shock were de- 
duced using ideal shock relations for a gas of y = 1. 4. For initial 
conditions of room temperature and atmospheric pressure with hydrogen con- 
centrations of 12 to 16 percent in air, the Von Neumann spike pressures 
and temperatures are in the range of 16 to 20 atmospheres and 1100° to 
1300° K. Delays have been computed from Eq. (30) with rate constants 
for the initiation reaction and reactions I, II, and III from refer- 
ence 4. Rate constants for reaction VI were obtained from Eq. (26) and 
the constants for reaction XI were taken as 

k -| i = 5.4x10® exp(-24, 000/RT) liters mole"-*-sec"^ 

which is simply 4.5 times the value suggested by Voevodsky.^® 

Results are shown in Fig. 4. Clearly there is an abrupt decrease 
in ignition lag around a hydrogen mole fraction of 0. 14; the ignition 
delay decreases by two orders of magnitude as the hydrogen concentration 
is increased 0.01. This provides further support for Belles' postulate 
that the condition 2kg > kgc^ must be satisfied for stable detonations 
in the hydrogen-oxygen system. 

Experimentally, the lean limit of detonability of hydrogen- air mix- 
tures lies in the range of 14 to 15 percent of hydrogen by volume . ^ 

CONCLUDING REMARKS 

It is hoped that the examples discussed in the preceding section il- 
lustrate the usefulness of these analytic solutions as applied to a 
variety of problems. Further analysis of this sort may help uncover more 
detailed information on the hydrogen- oxygen reaction. Some other applica- 
tions suggest themselves: 
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1. The concentrations of labile intermediates - H, OH, 0, HO2, etc. - 

ft 

can be computed as a function of time throughout the induction period. 

Thus, one can explore the feasibility of detecting these species experi- 
mentally - for instance, spectoscopically, or with EPR, or the mass 
spectrometer. 

2. By analyzing experimental ignition delay data obtained under a 
variety of conditions, it should be possible to glean information as to 
the mechanism of initiation, that is, the process by which atoms and 
free radicals are first generated from the reactants (by analysis of the 
term B in Eq. (24)). 

3. By computing delays near the low-pressure limits of detonability, 
it may be possible to gain insight as to the degree of coupling between 
shock front and reaction zone required in a stable detonation. 

4. If need be, other chemical reactions may be added to the scheme. 

(This may be desirable particularly in the high-pressure - low- temperature 
region. ) There is a restriction, however: reactions must be first order 

in the intermediates or products. Thus, reactions such as 

H + HgOg ►HgO + OH cannot be considered. 

5. Finally, the general method of attack used here should be useful 
in obtaining solutions to the kinetics for other isothermal branched chain 
explosions. 
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TABLE I. - NUMERICAL COMPARISON OF VARIOUS 


APPROXIMATIONS TO THE PARAMETER A 

[5 Percent Hg in air; rate constants 
from ref. 4. ] 


Steady- 

A '-l 

Temperature, °K 

state 
in - 

sec x 

1100 

1900 



Pressure, atm 



0.2307 

0. 4873 

OH 

i — i 
r< 

2.67X10 4 

3. 70x10 5 


*2 

2.39 

2.02 


(A 1 A 2) 1 / 2 

2. 526 

2.734 

H 

V 

10. 1x10 4 

4. 58x10 5 


A 2 

0. 6S5 

1. 616 


(AiAg) 1 / 2 

2. 630 

2. 720 

0 

1 — 1 
f< 

8.76X10 4 

16. 8x10 5 


• A 2 

0.793 

0.360 


( A 1 A 2) 1//2 

2. 638 

2.460 


True value 

2. 50 3X10 4 
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TABLE II. - COMPARISON OF NUMERICAL AND APPROXIMATE 
ANALYTIC CALCULATIONS OF HYDROGEN- 
AIR IGNITION DELAYS 


Temper- 

ature, 

°K 

Pressure, 

atm 

Mole frac- 
tion of H2 

Ignition delay, |isec 

Approximate 
(this paper) 

Numerical 

973 

1.0 

0.296 

240 

260 (ref. 1) 

1667 

1.0 

.296 

5.3 

5.1 (ref. 1) 

1100 

.2307 

.05 

866 

866 (ref. 4) 

1900 

.48 7 3 

.05 

33.6 

32.6 (ref. 4) 
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Fig. 2. - Effect of hydrogen concentration on hydrogen-air ignition 
delays^ at 594° C and 1 atmosphere. 




Fig. 3. - Effect of hydrogen concentration on hydrogen- 
air ignition delay 12 for 658 ±2° C and 1 atmosphere. 



Fig. 4. - Ignition delays for the Von Neumann spike condition 
in Chapman-Jouguet detonations of hydrogen -air in the 
neighborhood of the lean limit of detonability. 


